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Abstract 

Ah initio calculations show an antiferromagnetic-ferromagnetic phase transition 
around 9-10 GPa and a magnetic anomaly at 12 GPa in BiFeOs. The magnetic 
phase transition also involves a structural and insulator-metal transition. The G- 
type AFM configuration under pressure leads to an increase of the y component and 
decrease of the z component of the magnetization, which is caused by the splitting 
of the (1^2 orbital from doubly degenerate eg states. Our results agree with recent 
experimental results. 

Key words: G-type AFM structure;Phase transition;Dzyaloshinskii-Moriya 
interaction(DMI) 

PACS: 75.30.Et,75.30.GW,71.15.Mb 



1 Introduction 

Multiferroic materials have more than two of ferroelectric/antiferroelectric, 
ferromagnetic/ antiferromagnetic, and even ferroelastic/ antiferroelastic order- 
ing in same phase p!|2][3]H] . This feature makes them have potential applica- 
tions in information storage, spintronics,and sensors. The study on BiFeOa, a 
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multiferroic, which possess weak ferromagnetism and ferroelectric properties 
simultaneously at room temperature, has been recovered recently. [5|6|7] . It has 
long been known to be ferroelectric with a Curie temperature of about 1103 K 
and G-type antiferromagnetic(AFM) with a Neel temperature of 643 K. The 
large difference between the magnetic and ferroelectric ordering temperature 
makes the linear magnetoelectric coefficient smaller. Moreover,the -R3c space 
group permits a spiral spin structure in which the AFM axis rotates through 
the crystal with a long-wavelength period of 620A, which further reduces the 
observed G-type AFM magnetization. Fortunately, the decrease of magneti- 
zation can be compensated by doping in B-site of perovskites and fabricating 
thin film samples [5|9|TU] . The ferroelectricity is produced by the Bi-6s stereo- 
chemically active lone pair , which can only occur if the cation ionic site have 
broken inversion symmetry, while the weak magnetism is mainly attributed 
to Fe^+ ions. The coupling between magnetic and ferroelectric parameter is 
weak because they are driven by different ionic sites, and this agrees with 
the fact of large difference between the Curie temperature and AFM Neel 
temperature. Through first-principles calculations, we have shown that the 
rotation of the oxygen octahedra(antiferrodistortive(AFD) distortion) cou- 
ples with the weak ferromagnetism due to the Dzyaloshinskii-Moriya inter- 
action(DMI), considering the spin-orbital(SO) coupling effect and the non- 
coUinear spin configuration |11]. 

The study of BiFeOs under pressure should give us much insight in the 
understanding of the magnetoelectric coupling because it involves interest- 
ing structural and magnetic changes under pressure. Another metal-insulator 
transition and magnetic anomaly have already been reported recently [T2|T3] . 
Meanwhile, A. J. Hatt and coworkers performed the first principles calcula- 
tions on strain induced phase transition of BiFeOs film on the (OOl)-oriented 
substrate [HITS] . An isosymmetric phase transition accompanying with dra- 
matic structural change has been found. The two isosymmetric phases have 
the same space-group symmetry due to the constraints caused by coherence 
and epitaxy. The isosymmetric transition also indicates the coexistence of 
the rhombohedral and tetragonal phases. The rhombohedral film tends to be 
expanded with large compressive strain. And the lattice expansion can be 
accommodated by the increase of tetragonal domain. Further theoretical cal- 
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culations about the phase transition of BiFeOs need to be done to explain the 
detailed mechanism under pressure. We have performed first-principles calcu- 
lations to investigate the magnetic and electronic properties of BiFeOs under 
lower pressure range corresponding to the experiments to shed light on the 
mechanism. The remainder of this article is structured as follows: In section 2, 
the computational details of our calculations are given. We discuss the results 
in section 3. Our conclusion are given in section 4. 

2 Computational details 

We performe calculations within the local spin density approximation(LSDA) 
to DFT using the ABINIT package [T61IT7] . The ion-electron interaction is mod- 
eled by the projector augmented wave (PAW) potentials [T81IT9] with energy 
cutoff of 500 eV. We treat Bi 5d, 6s, and 6p electrons, Fe 4s, 4p,and 3d elec- 
trons, and O 2s and 2p electrons as valence states. 10 x 10 x 10 Monkhorst-Pack 
sampling of the Brillouin zone are used in calculations. LSDA-I-U method is 
introduced, where the strong coulomb repulsion between localized d states has 
been considered by adding a Hubbard-like term to the effective potential [20|21|22j . 

3 Results and discussion 

3. 1 FM behavior 

According to the experimental results and previous calculations P^P5] . we 
have constructed three structures, rhombohedral, monoclinic, and orthorhom- 
bic structures, with space group -R3c, Cm, and Pnma, respectively. The lattice 
parameters are given in Table 1 . The phase transition under pressure have been 
discussed elsewhere, and we suggest there exists a first-order phase transition 
from R3c to Cm structure around 9-10 GPa, and gradually to pure Pnma 
structure around 12 GPa which is close to the experimental values [T3] .We only 
report the magnetic and electronic behavior in this work. FM spin configura- 
tion for each structure are used to analyze the magnetization under pressure. 
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The FM magnetization are shown in Fig.l. One can see the spin value for the 
monochnic structure is generally larger than the other two structures, and it 
exhibits a decrease of spin magnetization with increasing pressure and leads to 
a flat curve after passing 9-10 GPa. The monoclinic and rhombohedral struc- 
tures both have the same spin value at all pressure range, indicating they have 
same structure dependent spin configuration. Moreover the phase transition 
between these two can be much easier than others. The transition between 
these two structures happens at 9-10 GPa. That between rhombohedral and 
orthorhombic structures take places at higher pressure of 12 GPa. 



3.2 AFM-FM spin order transition 



Meanwhile, for rhombohedral structure we set up three spin configura- 
tions, FM, AFM, and G-type AFM spin structures. The total energy per unit 
cell for different spin configuration have been calculated and shown in Fig. 2, 
as well as considering the on-site Coulomb interaction term, U. In LSDA-I-U 
calculations, U and J are defined as 

U = — 777 V < m,m'\Vee\m,m' >= F^, (1) 



1 ^2 _^ ^4 

2/(2/+ 1) , ' I eel , , y ) 

\ ' m^m' ,m' 

where Vee are the screened Coulomb interaction among the nl electrons. F^, 
F"^, andF^ are the radial Slater integrals for d electrons in Fe. It is apparent 
that there are intersection points between the FM and AFM energy curves. As 
pressure approaches the intersection points, both LSDA+U and LSDA give 
lower AFM energy. However FM energy are much favorable after pressure ex- 
ceeds the intersections, indicating an AFM-FM phase transition around the 
pressure value of 9-10 GPa. The LSDA+U can produce relatively higher en- 
ergy. On-site coulomb interaction , U is the energy needed to put two electrons 
in the same site. The value of U among the electrons in transition-metal d or- 
bitals are one magnitude higher than the Stoner parameter. An appropriate 
band gap can be obtained in transition-metal oxides with properly choosing 
the U parameter . In our calculation we choose U=2 eV as it is rightly un- 
der the critical value to preclude the DMI in G-type AFM spin configuration. 
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where the Fe ions are arranged antiferromagnetically along x direction |24]. 
We take into account the non-colhnear spin structure and spin-orbital(SO) 
interaction in G-type configuration. One can see that the G-type spin struc- 
ture leads to a lower energy comparing with the FM and AFM one within 
LSDA+U scheme. Except the anomaly around the first transition pressure of 
9-10 GPa, there is another one around 12 GPa. The relaxation results show 
the structure is changed into an orthorhombic phase at this pressure. 

3.3 Exchange interaction 

Consideration the Heisenberg model, 

A£;=-l/2 5]J,„,Si-Sj, (3) 

for FM spin configuration, the total energy involving the spin exchange 
interaction can be written as, 

EpM = Et — ZcJzS'^, (4) 

where Et is the total energy without the spin, Zc is the number of nearest 
neighboring Fe ions, and is the exchange integral. While the AFM total 
energy has the form, 

Eafm = Et + ZcJzS"^. (5) 

Therefore, we can determine the exchange parameters from the energy differ- 
ence between different spin configurations. We set up the FM structure with 
spin direction along z axis while the G-type along x direction. From Fig. 3, 
it is obvious that R?>c space group tends to produce a much favorable G- 
type structure. LSDA+U gives lower AFM exchange integral than LSDA. It 
is worth pointing out again that the anomaly takes place around the criti- 
cal pressure value. Exchange integral even lies in the same level in these two 
schemes. G-type structure has lower exchange interaction under the whole 
pressure range. Two anomalies can be found at 9 and 12 GPa, respectively. 
It is consistent with energy calculations. When pressure exceeds the critical 
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value, AFM exchange integral is changing into a positive value where it ap- 
plies to FM spin configuration. This shows the AFM-FM transition occurs at 
the critical pressure accompanying with the structural phase transition which 
agrees well with the previous total energy calculations. The phase graph under 
pressure is shown in Fig. 4. We suggest the rhombohedral structure maintains 
before the critical pressure value of 9 GPa. A combination of these three 
structures exist between 9 and 12 GPa, accompanying with an AFM-FM spin 
transition. A pure orthorhombic phase will be found after 12 GPa while the 
FM spin structure remains. 



3.4 G-type AFM vectors 



In G-type spin structure, we take into account the non-coUinear and SO 
coupling effect. The AFM vectors under pressure are reported in Fig. 5. It is 
apparent that the AFM spin in x direction cancel out and a resultant mag- 
netization along y and z direction can be obtained due to the DMI. DMI is 
caused by the interaction of neighboring Fe sites which can be described by 



-^Fil,Fe2 - JFel,Fe2(Sl ' §2) + D^el,Fe2(Sl X S2) + S{R) ■ T'§li^Fe2 ' ^2, (6) 

in the second order perturbation calculation . The first term on the right 
hand side of the Eq. (6) corresponds to the usual isotropic superexchange 
interaction, and the second term is the DMI. The Hamiltonian for the system 
reads. 



HsiFeOs ~ ~2 ^ Jli,2iSij • S2j + ^ Dij^2jSij X 82^. (7) 

<li,2j> <li,2j> 

The first term is the symmetric superexchange, and the second one is the 
antisymmetric DMI contribution. Jit,2j is a constant similar to the exchange 
interaction. D is the DMI constant and determined by the sense of rotation of 
the neighboring oxygen octahedra. D reads by the second order perturbation 
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in the case of one electron per ion 

B%,Fe2 = i^^/U)[bnniFel-Fe2)Cn'niFe2-Fel)-Cr,niFel-Fe2)br,'n{Fe2-Fel)^^ 

(8) 

where U is the energy required to transfer one electron from one site to its 
nearest neighbor, a parameter similar to on-site Coulomb interaction in our Ab 
initio computation, and inversely proportional to D. The spin value along z 
direction is depressed after pressure exceeds the critical value while spin mag- 
netization along y direction has the opposite trend and increases with pressure 
increasing. Firstly, the net magnetization has components along z and y direc- 
tion simultaneously and deviate away from the z direction as pressure exceeds 
12 GPa, resulting in zero component along z, while that along y direction 
increases and maitains a constant value after 12 GPa. The magnetization per 
unit cell is calculated based on the LSDA+U method and it is underestimated 
in this way. Therefore greater value is expected in LSDA calculations [21]. 

3. 5 Electronic properties 

In order to shed light on the electronic properties under pressure unam- 
biguously, the total density of states(DOS) before and after exerting pressure 
are given in Fig. 6. The orbital resolved DOS(ODOS) for Fe 3d orbitals are 
given in Fig. 7 and Fig. 8, respectively. From Fig. 6, it can be seen that a 
semiconducting band gap of 1.7 eV is produced within LSDA+U while the 
band gap vanishes under pressure of 12 GPa, suggesting an obvious insulator- 
metal (IM) transition. The IM transition is mainly caused by the shift of the 
states of Fe 3d electrons in the vicinity of fermi energy. The finite DOS of 
these electrons cut through the fermi level and form a FM spin structure un- 
der 12 GPa. From Fig. 7, one can see all up-spin electrons are occupied and 
down-spin electrons are in conduction band. While almost all up and down 
spin states are partially filled under pressure in Fig. 8. Three- fold degenerate 
states t2g coming from dxy, dyz, and d^z orbitals remain degenerate, while d^^ 
orbital splits from the two- fold degenerate states Cg. It is this orbital that 
makes Cg-eg AFM interaction reduced. The splitting of the orbital under pres- 
sure is significant, suggesting the complete depression of AFM interaction and 
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the occurrence of FM spin structure. This is also the reason for the decrease 
of magnetization of G-type AFM structure along z direction under pressure. 
Meanwhile the down-spin Fe Sc? states under pressure cut across the fermi 
energy and lead to the conducting behavior. 

4 Conclusion 

The total energy, magnetic and electronic properties of BiFeOs under 
pressure are calculated based on the LSDA and LSDA+U scheme. Results 
show two anomalies can be found at 9-10 GPa and 12 GPa, respectively. The 
first one is the critical pressure for first-order phase transition accompanying 
with AFM-FM transition. Meanwhile the behavior under critical pressure also 
involves an IM transition. The second one does not involve further magnetic 
transition but structural transitions. The magnetization of the G-type AFM 
spin structure in y direction increases while the z direction component de- 
creases, which can be explained by the splitting of the (i^2 orbital from doubly 
degenerate tg states. 
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Table 1 

Current used lattice parameters for rhombohedral, monoclinic, and orthorhombic 
structures 





R3c 


Cm 


Pnma 


a(A) 


5.459 


5.7900 


5.5849 


b(A) 




5.6899 


7.6597 


c(A) 




4.1739 


5.3497 


an 


60.36 










91.99 




v(A3) 


115.98 


137.42 


113.12 


Bi 


(2a):0,0,0 


(2a):0.9376,0,0.0685 


(4c):0.0536,0.2500,0.9886 


Fe 


(2a):0.2308,0.2308,0.2308 


(2a):0.5110,0,0.4961 


(4b):0,0,0.5 





(6b):0.5423,0.9428,0.3980 


(2a):0.5626,0,0,9489 


(4c):0.9750,0.2500,0.4060; 






(4b):0.7958,0.7603,0.4231 


(8d):0.2000,0.9540,0.1945 
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Figure captions: 

Fig.l FM magnetization for tliree pliases of BiFeOs under pressure. 
Fig. 2 Tlie total energy as functions of pressures. 

Fig. 3 The exchange integrals for different spin structures as functions of 
pressure. 

Fig. 4 The phase graph of BiFeOa under pressure. 

Fig. 5 G-type AFM vectors variations with respect to pressure. 

Fig. 6 Total DOS under ambient and transition pressure. 

Fig. 7 ODOS for Fe dxy,dyz,dz2,dxz, and d^^-yi orbitals under ambient 
pressure. 

Fig. 8 ODOS for Fe dxy,dyz,dz2,dxz, and dx2-y2 orbitals under transition 
pressure. 
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Fig. 1. FM magnetization for three phases of BiFeOs under pressure. 
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Fig. 2. The total energy as functions of pressures. 
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Fig. 3. The exchange integrals for different spin structures as functions of pressure. 
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Fig. 4. The phase graph of BiFeOa under pressure. 
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Fig. 5. G-type AFM vectors variations with respect to pressure. 




Fig. 6. Total DOS under ambient and transition pressure. 
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